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Abstract 

In this paper we are interested in the mathematical and 
numerical analysis of the time-dependent Galbrun equa- 
tion in a rigid duct. This equation models the acoustic 
propagation in the presence of a flow 1 1 1. We propose a 
regularized variational formulation of the problem, in the 
subsonic case, suitable for an approximation by Lagrange 
finite elements, and corresponding absorbing boundary 
conditions. 

1 Introduction 

, Sound propagation in a flow is a subject of great inter- 
est for numerical analysis. The main application concerns 
noise reduction in aeronautics. The understanding of phe- 
nomena of interaction between acoustic waves and flows 
is a crucial feature to find the components that efficiently 
reduce the sound. 

Most of the mathematical and numerical studies for 

' transient linearized aeroacoustics are based on linearized 
Euler equations, whose unknowns are the perturbations 

■ of velocity and pressure. Nevertheless, Galbrun equa- 

', tions, whose unknown is the lagrangian perturbation of 
displacement, are an attractive alternative to model the 
phenomenon of acoustic propagation in a flow. Indeed, 

, these equations seem to have a structure similar to equa- 
tions in electromagnetism and elastodynamics and they 
allow a simple treatment of boundary conditions. 

Recent studies on Galbrun equations have focused on 

! time harmonic dependence. We are interested here in 
studying the transient case. 

This problem raises many theoretical and numerical 

; difficulties. The major one is the lack of a natural func- 
tional frame for a variational setting of the problem. 
Moreover, a naive numerical resolution of the problem us- 
ing standard nodal finite elements for the space discretiza- 
tion is unstable (see Figure©. In order to overcome these 
difficulties , we propose a regularization method analo- 
gous to the one developed for Maxwell's equations 
This allows to derive a variational formulation suitable 
for an approximation by Lagrange finite elements. 

2 Statement of the problem 

Let = Mx] — /i,/i[bea two-dimensional rigid duct con- 
taining a compressible fluid. We suppose that this fluid is 



in motion and that the flow can be described by a profile 
of Mach's number M{y) {M g] - 1, 1[). If / denotes 
the source of acoustic waves, we have to solve the tran- 
sient Galbrun equation in dimensionless form ( ^ is the 
perturbation of lagrangian displacement): 
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where gj 



|j + M{y)-^ is the material derivative and 
n the unit outward normal vector to dO,. We complete the 
equations Q-© by adding initial conditions. 



3 The Lagrangian vorticity equation 

When / e i7(curl, Q,), applying the curl operator to (0, 
we prove that the curl^ = ip satisfies: 

dI = ^^^'^^^D^(^) +™'^-^ innxRX (3) 

Note that when the flow is uniform {M{y) is constant), 
the vorticity ijj can be computed independently of ^ and 
the solution is given by : 



V'lx, y, t) = a{x - Mt, y) + xp{x - Mt, y) + 
~ a)(™rl/)(a, y, t - ^)da 



where a et /3 are two functions that depend only of initial 
conditions of the problem. 

4 Regularized Formulation Galbrun equation 

The idea of the regularization, initially introduced for 
Maxwell's equations, was extended to the time har- 
monic Galbrun equation by A. S Bonnet-Bendhia and al 
(2001). The idea consists in adding the artificial term 
s curl (curl —^fj) to Galbrun equation. 
We replace the initial value problem Q-© by the equiv- 
alent regularized problem : 

§^ - V(div|) + scurl(curl$-'(/') = / in x R+ 



m (^-2Ar(y)^i^) =curl/ in fix 



g,(^-2Ar(y)f 
^ • n = 0, curl^ = ijj 

where s is a non negative parameter. 
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In this paper we restrict to the case of uniform flow 
{M'{y) =0 Vy s] — /i, h\) for which the problem in ip 
and are decoupled. The study of the coupled system 
will be a subject of future works 

5 Mathematical analysis of the regularized problem 

When the Mach number is constant the vorticity ijj is 
known and the problem in ^ is written : 



V(div^) + scurl (curl^) 
^ ■ n = 0, curl^ = ip 
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where fs = f + s curl ip. 

The second boundary condition (curl^ = ip) of system 
^ is necessary for the equivalence with the initial prob- 
lem Q-©. The property of ellipticity of the spatial op- 
erator — V(div) + scurl (curl) (which equals —A if 
s = I) allows to carry out the mathematical and numeri- 
cal study of this problem in a classical frame. 
We consider H := Ho(l^) x ^2(17)2 where 

Uo{n) := e {H\n)f / ^■n = 0,ondn}. 

To apply Hille-Yosida theorem, we introduce a new un- 
known C = D^/Di- If we pose U = (^, then we 
can rewrite ^ under the following form : 




dx 



(6) 



y-V(div|) + s curl (curl|) + M-^j 

The domain of the unbounded operator As is defined by : 

D{As) = {U={tCY^n tel que AsU G . 

Using Costabel's identity |3| and Hille-Yosida's theory, 
we prove that the operator Ag is maximal monotone 

Theorem 5.1 .• // min(l, s) > M^, then for fs G 
C^(IR+; ^2(17)2) and sufficiently regular initial data, 
problem (ISj has a unique solution which satisfies: 

i € Ci(M+; Ho) n C2(M+; l2(17)2) 



6 The absorbing boundary conditions 

For solving numerically the problem, we need to truncate 
the unbounded domain Q.. The truncated domain ilf, := 
] — R; i?[x[— /i; h[ requires the introduction of absorbing 
boundary conditions (ABCs) on the artificial boundaries 
r± := {(x,y) e 17, X = ±R} {R > 0). The difficulty 
is to find appropriate boundary conditions adapted for the 
regularized formulation. For s = 1, the boundary condi- 
tions that we propose are the following : 
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where n = (±1, 0) is the unit outward normal vector of 
and r = (0, ±1). The main properties of conditions 
Q are 

i) These are exact conditions for y-independent solu- 
tions (plane waves) : this is why we speak of first 
order conditions. 

ii) Well-posedness: the truncated problem is well 
posed and for / = we have the energy decay re- 
suh: 

d 
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iii) These are compatible with a variational formula- 
tion of the regularized problem, namely: 



with (•, ■) is the L2-inner product and 
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Remark 6.1 The construction of good absorbing bound- 
ary conditions is not trivial. For example the following 
natural ABC's : 



= 



r± 



r± 



(9) 
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are still exact for plane waves and variational, but not 
stable. 

7 Discretization of Galbrun equations 

The Lagrange finite element approximation of ^ leads 
to tlie following ordinary differential system : 



+ 



where M.h is the mass matrix, and A/j, B/j, C^"^ and Bj^* 
are the matrices respectively associated to the bilinear 
forms a(., .), b{., .), c^^{., .) and (f^{., .). 

For the time discretization of (fTTT) we use a centered 
second order finite difference scheme : 

/■n+l od-n I tn— 1 ^n+l 
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8 Numerical Simulation 

In this first experiment, we simulate a wave initially ex- 
cited by a rotational source located in the center of the 
domain, in the presence of a horizontal uniform flow with 
M = 0.5. The first simulation shows that the method is 
not stable if the equation is not regularized {s = 0). 






Figure 1 : Euclidean norm of ^ at at ti 
1.75s and ts = 2s. Case s = 



1.5s and t2 



The second simulation corresponds to the regularized 
case {s = 1). When can distinguish the two parts of 
the Lagrangian displacement : the irrotational part corre- 
sponds to the outer circular wavefront (whose radius in- 
creases with time) while the rotational part corresponds to 
the inner circular wavefront (whose radius remains con- 
stant). Both are centered at a point which is convected by 
the flow. 




Figure 2: Euclidean norm of ^ at ti 
1.75s and t^, = 2s. Case s = 1 



1.5s and t2 



In the second experiment, we simulate Galbrun wave 
propagation in an infinite rigid duct. The wave is excited 
by a Gaussian signal in time and a quasi-punctual irrota- 
tionnal source in space in the uniform flow with M = 0.5. 
We remark on this result that the ABC prevents (partially) 
the unphysical reflexion. 




Figure 3: Evolution of the EucUdean norm of the dis- 
placement ^ 
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